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ABSTRACT 

We construct size distributions for carbonaceous and silicate grain populations in dif- 
ferent regions of the Milky Way, LMC, and SMC. The size distributions include sufficient 
very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) 
to account for the observed infrared and microwave emission from the diffuse interstellar 
medium. Our distributions reproduce the observed extinction of starlight, which varies 
depending upon the interstellar environment through which the light travels. As shown 
by Cardelli, Clayton & Mathis in 1989, these variations can be roughly parameterized 
by the ratio of visual extinction to reddening, Ry ■ We adopt a fairly simple functional 
form for the size distribution, characterized by several parameters. We tabulate these 
parameters for various combinations of values for Ry and be , the C abundance in very 
small grains. We also find size distributions for the line of sight to HD 210121, and for 
sightlines in the LMC and SMC. For several size distributions, we evaluate the albedo 
and scattering asymmetry parameter, and present model extinction curves extending 
beyond the Lyman limit. 

Subject headings: dust — extinction — ISM: clouds — Magellanic Clouds 



1. Introduction 

Mathis, Rumpl, & Nordsieck (1977, MRN) constructed their classic interstellar dust model on 
the basis of the observed extinction of starlight for lines of sight passing through diffuse clouds. 
Strong absorption is observed at 9.7 and 18^m, corresponding to stretching and bending modes 
in silicates. The strong extinction feature at 2175 A can be approximately reproduced by small 
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graphite particles (Stecher Sz Donn 1965; Wickramasinghe & Guillaume 1965). The simplest model 

incorporating both silicate and graphite material consists of two separate grain populations, one of 
silicate composition and one of graphite composition. MRN found that the extinction curve (i.e. 
the functional dependence of the extinction on the wavelength A) is well reproduced if the grain 
size distribution (with identical form for each component) is given by 

drigr = Cn]ia~^-^da, Omin < a < Omax (1) 

with amin = 50A and Omax = 0.25/xm; ngr(a) is the number density of grains with size < a and nn 
is the number density of H nuclei (in both atoms and molecules). MRN adopted spherical grains, 
for which Mie theory can be used to compute extinction cross sections, and we shall do the same; 
in this case a is the grain radius. Drainc & Lee (1984) extended the wavelength coverage of the 
MRN model, constructed dielectric functions for "astronomical silicate" and graphite, and found 
the following normalizations for the size distribution: C = 10"^^-^^ (10~^^'^^) cm^-^ for graphite 
(silicate) . 

Since the development of the MRN model, more observational evidence has become available; 
some of these new observations require revisions of the model. First, the extinction curve has been 
found to vary, depending upon the interstellar environment through which the starlight passes. 

Cardelli, Clayton, & Mathis (1989, CCM) found that this dependence can be characterized fairly 
well by a single parameter, which they took to be Ry = A{V)/E{B — V), the ratio of visual 
extinction to reddening. CCM have fitted the average extinction curve A{\)/A{V) as functions 
of A and Ry- For the diffuse ISM, Ry ~ 3.1; higher values are observed for dense clouds. Kim, 
Martin, &; Hendry (1994) used the maximum entropy method to find smooth size distributions, 
for silicate and graphite grains, for which the extinction for Ry = 3.1 and 5.3 is well reproduced. 
Their Ry = 5.3 distribution has significantly fewer "small" grains (a < 0.1/im) than their Ry = 3.1 
distribution, as well as a modest increase at larger sizes. This result was expected, since generally 
there is relatively less extinction at short wavelengths (provided by small grains) for larger values 
of Ry. 

Observations of thermal emission from dust have provided another challenge to the MRN 
model. Emission in the 3 to 60/xm range, presumably generated by grains small enough to reach 
temperatures of 30 to 600 K or more upon the absorption of a single starlight photon (see e.g. 
Draine & Anderson 1985), imply a population of very small grains (with a < 50 A). The non- 
detection of the 10/um silicate feature in emission from diffuse clouds (Mattila ct al. 1996; Onaka et 
al. 1996) appears to rule out silicate grains as a major component of the a < 15 A population (but 
see note added in proof). Emission features at 3.3, 6.2, 7.7, 8.6, and 11.3//m (see Sellgren 1994 for a 
review) have been identified as C-H and C-C stretching and bending modes in polycyclic aromatic 
hydrocarbons (Leger & Puget 1984), suggesting that the carbonaceous grain population extends 
down into the molecular regime. Recent observations of dust-corrclated microwave emission has 
been attributed to the very small grain population (Draine &; Lazarian 1998a). 

The abundance of very small grains required to generate the observed IR emission from the 
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diffuse ISM is not yet well-known. In the model of Desert, Boulanger, & Puget (1990), polycyclic 
aromatic hydrocarbon (PAH) molecules with less than 540 C atoms (equal to the number of C 
atoms in a spherical graphite grain with a ~ 10 A) lock up a C abundance^ of ~ 4 x 10^^. Li & 
Draine (2001) compare observations of diffuse galactic emission with detailed model calculations 
for grains heated by galactic starlight and find that a C abundance ~ 4-6 x 10~^ is required in 
hydrocarbon molecules with < 10^ C atoms. They conclude that the emission is best reproduced 
if the very small grain population is the sum of two log-normal size distributions:^ 
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l + erf[3o-/\/2 + ln(ao,i/3.5A)/a\/2] ' 

where mc is the mass of a C atom, p = 2.24 gcm~^ is the density of graphite, 6c, i = 0.755c, ^c,2 = 
0.256c 5 6c is the total C abundance (per H nucleus) in the log-normal populations, ao,i = 3.5 A, 
00,2 = 30 A, and a = 0.4. 

Draine &: Lazarian (1998b) estimated the electric dipole radiation from spinning grains, and 
found that it could account for the dust-correlated component of the diffuse Galactic microwave 
emission if 6c ~ 2 x 10~^. More recent modelling confirms that the microwave emission can be 
reproduced with 6c ~ 2-4 x 10"^ (B. T. Draine 8z A. Li 2001, in preparation). 

Our goal here is to find size distributions which include very small carbonaceous grains^ (in 
numbers sufficient to explain the observed infrared and microwave emission attributed to this 
population) and are consistent with the observed extinction, for different values of Ry in the 
local Milky Way and for regions in the Large and Small Magellanic Clouds. We consider several 
values of 6c, since the C abundance in very small grains is not yet established. We discuss the 
observational constraints and our method for fitting the extinction in §2, present results in §3, and 
give a discussion in §4. The size distributions obtained here will be employed in separate studies, 
including an investigation of photoelectric heating by interstellar dust (Weingartner Sz Draine 2001). 



^By "abundance", we mean the number of atoms of an element per interstellar H nucleus. 

^The log-normal distribution with ao,i ~ 3.5 A is required to reproduce the observed 3-25/im emission, and the 
= 30A component is needed to contribute emission in the DIRBE 60/iin band. 

^We take "carbonaceous grains" to refer to graphitic grains and PAH molecules. Although not all carbonaceous 
grains arc graphite, we will continue to refer to the dust model considered here as the "graphite/silicate" model, for 
simplicity. 
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2. Fitting the Extinction 

2.1. "Observed" Extinction 

For the "observed" extinction Aa]~,s{Ry, A), we adopt the parametrization given by Fitzpatrick 
(1999). Bohhn, Savage, &: Drake (1978) found that the ratio of the total neutral hydrogen column 
density A^h (including both atomic and molecular forms) to E{B — V) is fairly constant for the 
diffuse ISM, with value 5.8 x 10^^ cm~^. This provides the normalization for the extinction curve: 
A{V)/N}i = 5.3 X 10"^^ cm^. The normalization is less clear for dense clouds, because of the 
difficulty in measuring iVn.^ CCM found that A{X)/A{I) appears to be independent of Ry for 
A > 0.9^m (= / band), suggesting that the diffuse cloud value of A{I)/N}i = 2.6 x 10~^^ cm^ may 
also hold for dense clouds (see, e.g., Draine 1989); we adopt this normalization. 



2.2. Functional Form for the Size Distribution 

Lacking a satisfactory theory for the size distribution of interstellar dust, we employ functional 
forms for the distribution which (1) allow for a smooth cutoff for size a > at, with control of the 
steepness of this cutoff; and (2) allow for a change in the slope dlnngr/dln a for a < at- We adopt 
the following form: 

1 dugi _ -p./ \ f a . , J 1 , 3.5A < a < at,g 



a \at^gj [exp|-[(a-at,g)/a 



nn da a V«t,g/ ' I exp {-[(a - at,g)/ac,g]'^| , a > at,g 

(4) 

for carbonaceous dust [with D{a) from eq. (2)] and 

1 drigr f a ' 3.5A<a<at,s 



nn da a V«t,s/ ' ^' I exp {-[(a - at,s)/cic,s]^} , a > at,s 

for silicate dust. The term 

[(l-Z^/a/ot) , /? < 

provides curvature. The form of the exponential cutoff was suggested by Greenberg (1978). The 
structure of the size distribution -D(a) for the very small carbonaceous grains has only a mild effect 
on the extinction for the wavelengths of interest; we adopt the same values as Li & Draine (2001) 
for ao,i = 3.5 A, ao,2 = 30 A, and a = 0.4, and the same relative populations in the two log-normal 
components (6c, 1 = 0.756c) bc,2 = 0.256c)) but will consider different values of 6c- Thus equation 
(4) has a total of six adjustable parameters (6c, Cg, at,g, ac,g, ag, /3g), with another five parameters 
(Cs, o,t,fi, Qc.s- fis^ fts) ill equation (5) for the silicate size distribution. 



''Kim & Martin (1996) compiled a set of sight lines for which both A{V)/Nyi and Rv are observationally deter- 
mined. Their data are consistent with A(V)/Nh being independent of Rv, but the uncertainties are large. 
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2.3. Calculating the Extinction from the Model 

The extinction at wavelength A is given by 

A(A) = (2.57rloge)y"dlna^^^^^a3Qext(a,A) , (7) 

where A'gi.(a) is the column density of grains with size < a and Qcxt is the extinction efficiency 
factor, which wc evaluate (assuming spherical grains) using a Mie theory code derived from BHMIE 
(Bohren & Huffman 1983). 

We adopt silicate dielectric functions based on the "astronomical silicate" functions given by 
Draine &: Lee (1984) and Laor & Draine (1993), but differing in the ultraviolet. The "astronomical 
silicate" dielectric function e = e\+ie2 of Draine & Lee (1984), based on laboratory measurements of 
crystalline olivine in the ultraviolet (Huffman &: Stapp 1973), contains a feature at 6.5/xm~^. Kim & 
Martin (1995) have pointed out that this feature, which is of crystalline origin, is not present in the 
observed interstellar extinction or polarization. Wc have therefore excised this feature from £2 and 
"redistributed" the oscillator strength over frequencies between 8 and 10/um^^; we then recomputed 
ei using the Kramers-Kronig relation (Draine &: Lee 1984). (The resulting "smoothed astronomical 
silicate" dielectric functions are available at http://www.astro.princeton.edu/~draine/.) 

For carbonaceous grains, we adopt the description given by Li & Draine (2001), in which the 
smallest grains are PAH molecules, the largest grains consist of graphite, and grains of intermediate 
size have optical properties intermediate between those of PAHs and graphite. For PAHs, Li & 
Draine estimate absorption cross sections per C atom, for both neutral and ionized molecules. Li 
& Draine estimate PAH absorption near 2175 A by assuming that the 2175 A absorption profile is 
in large part due to the PAH population; our adopted PAH absorption cross sections near 2175 A 
therefore agree - by construction - with the observed 2175 A profile. We convert to a size-based 
description by assuming a C density p = 2.24gcm~^, and we assume that 50% are neutral and 
50% are ionized (the ionization state affects the absorption by these grains at A > 0.6/xm). We 
take graphite dielectric functions from Draine &: Lee (1984) and Laor & Draine (1993) and adopt 
the usual "1/3 — 2/3" approximation: Qcxt = [Qcxt(e||) + 2Qcxt(e±)]/3, where ey and e± are the 
components of the graphite dielectric tensor for the electric field parallel and perpendicular to 
the c-axis, respectively. Draine & Malhotra (1993) showed that the 1/3 — 2/3 approximation is 
sufficiently accurate for extinction curve modelling. 



2.4. Abundance/Depletion Constraints 

Given estimates of the abundances and interstellar depletions of the elements incorporated in 
dust and the mass densities of the grain materials, we can estimate the total volume per H atom, 
T4ot; in the carbonaceous and silicate grain populations. For a long time, solar abundances were 
used for this purpose (see Grevesse k. Sauval 1998 for a recent compendium of solar abundances). 
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Recent evidence, e.g. from measurements of abundances in the atmospheres of B stars, suggest that 
the abundances in the present-day ISM may be substantially lower than the solar values (sec Snow 
& Witt 1996; Mathis 1996, 2000; and Snow 2000 for reviews). However, Fitzpatrick & Spitzer 
(1996) concluded that S has solar abundance in the ISM, and Howk, Savage, & Fabian (1999) 
found solar abundances of Zn, P, and S along the line of sight to ^ Columbae. Thus, interstellar 
abundances are not yet well-known. 

We adopt the solar C abundance of 3.3 x 10~^ (Grevesse & Sauval 1998) and assume that 
30% is in the gas phase.^ With the ideal graphite density of 2.24 gcm~^, we find Vtot,g ~ 
2.07 X 10~^^ cm^ for carbonaceous dust. To estimate the total volume in amorphous silicates, 
we assume a stoichiometry approximating MgFeSiO.i, with mass number per structural unit of 172. 
Since Si, Mg, and Fe have similar abundances in the Sun and arc all highly depleted in the ISM 
(Savage & Sembach 1996), we simply assume that the Si abundance in silicate dust is equal to its 
solar value of 3.63 x 10~^. We adopt a density of 3.5gcm~^, intermediate between the values for 
crystalline forsterite (Mg2Si04, 3.21 gcm~^) and fayalite (Fe2Si04, 4.39 gcm~^). Thus, we estimate 
Vtot,s ~ 2.98 X 10-2^ cm^H-^ for silicate dust. 



2.5. Method of Solution 

For a given pair of values {Rv, be), we seek the best fit to the extinction by varying the powers 
Q!g and as', the "curvature" parameters /3g and (3s; the transition sizes at.g and at,s; the upper 
cutoff parameters Oc^g and Ocs! and the total volume per H in both the carbonaceous and silicate 
distributions, V^ot.g and Vtot.s, respectively. 

We use the Levenberg-Marquardt method, as implemented in Press et al. (1992), to fit the con- 
tinuous extinction between 0.35/[im~^ and 8//m~^.® We evaluate the extinction at 100 wavelengths 
Aj, equally spaced in In A, and minimize one of two error functions. In the first case (hereafter "case 
A") we minimize = Xi + Xy- 

The first term in gives the error in the extinction fit: 

2 (In^obs - In^mod)^ /o^ 

xi = 2^ ^ ' y^) 

i * 

where ^obs('^i) is the average "observed" extinction (§2.1), ^mod(-^i) is the extinction computed for 
the model [equation (7)], and the (jj are weights. When evaluating Amodi we verify that the integral 



^Cardelli et al. (1996) and Sofia et al. (1997) found a gas-phase C abundance of 1.4 x 10 larger than the 
w 1 X 10~* that we assume. 

®The lower limit of 0.35/[im~^ was chosen so as to avoid infrared absorption features, most notably the SAfjja C-H 
stretch feature. Extinction data for > 8/um~^ are very limited. 
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in equation (7) is evaluated accurately. We take the weights CTj =1 for l.l^um -"^ < A ^ < 8/xm ^ 
and CTj"^ = 1/3 for < l.l/xm~^, since the actual IR extinction is uncertain. 

The term Xy ^ penalty which keeps the total volumes in the carbonaceous and silicate grain 
populations from grossly exceeding the abundance/depletion-limited values found in §2.4. We take 

Xv = 0.4[max(yg, 1) - l]^-^ + 0.4[max(V;, 1) - l]^'^ , (9) 

where % = Viot,g/2.07 x IQ-^^cm^H-i and K = Vtot,s/2.98 x IQ-^^cm^H-^ 

Given our assumption that A{I)/Ny{ is independent of Ry, the extinction for higher Ry can be 
fit using less total grain volume. It seems highly unlikely that material is transferred from grains to 
the gas phase as gas and dust cycles into regions of higher density. Thus, we also consider a second 
case ("case B") for which the grain volumes are held fixed at approximately the values found for 
Rv = 3.1: Vtot.s = 3.9 x 10-^7 cm^H-^ and Vtot,g = 2.3 x IQ-^^cm^H-^ In this case, we seek to 
minimize Xi- 

3. Results 
3.1. Dust in the Milky Way 

3.1.1. Size Distributions and Extinction Fits 

We have generally found, in fitting the extinction, that x^ varies only slightly with the silicate 
cutoff parameter ac,s until ac,s exceeds a critical value of ^ 0.1/um (for Milky Way dust; see Figure 
1). As ac,s increases, the silicate grains contribute less short-wavelength extinction, and a large 
abundance of small carbonaceous grains is required to pick up the slack. When ac,s ^ O.l/im, the 
2175 A hump is overproduced. Although x^ is nearly constant for ac,s < 0.1/xm, it does increase 
slightly with ac,s- Consequently, our fitting algorithm returns very small values for ac,s, for which 
the silicate size distribution drops off very sharply at the large-size end. Since such sharp cutoffs 
are unlikely to occur in nature, we have opted to fix ac,s = 0.1/im. 

In Table 1 we list the values of the distribution parameters for which the extinction with 
Rv = 3.1, 4.0, and 5.5 is best fit, for various values of bc-^ These distributions are displayed in 
Figures 2 through 6. 

In Figure 7 we display ^obs and Amod for case A, the three values of Ry, and the highest values 
of 6c included in Table 1, in a log- log plot, to give a sense for the fit quality over the entire range of 
A~^. In Figures 8 through 12, we display extinction curves for be = and for the highest value of 
be included in Table 1; we show the contribution from each of the grain distribution components. 



These parameters [and a FORTRAN subroutine that returns dngr/da(o)] are also available in electronic form on 
the World Wide Web at www.cita.utoronto.ca/~weingart. 
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In Table 1 we also display x^, xf , and xl = Z^jCln^obs ~ In^mod)^- For a given value of Ry, 
the error functions do not vary substantially with be until a critical value of be is reached, at which 
point the error functions increase dramatically (see Figure 13). Clearly, extinction evidence alone 
does not constrain be well except that 5c < 6 x 10~^ for the Ry = 3.1 extinction law, 6c ^ 4 x 10~^ 
for Rv = 4, and fee ^ 3 x 10~^ for Ry = 5.5. In each case, the upper limit on be is reached when 
the very small carbonaceous particles account for 100% of the 2175A extinction feature. 

In assessing the quality of the extinction fits, one must bear in mind that (1) the dielectric 
functions used are certainly not correct in detail, even for bulk material, (2) the surface monolayers 
of grains are likely to differ from bulk materials, (3) the true size distributions undoubtedly differ 
from the adopted functional form, and (4) the interstellar grains are appreciably non-spherical. 
Therefore, a precise fit is not to be expected. One should also remember that the adopted PAH 
absorption cross section in the vacuum ultraviolet was constructed to fit the interstellar 2175 A 
profile, and the silicate dielectric function in the vacuum ultraviolet was modified to suppress 
structure not present in the observed interstellar extinction. 

3.1.2. Further Results 

Although neutral H gas is opaque for wavelengths shortward of the Lyman limit, extinction by 
dust at such wavelengths could have important observational consequences within ionized regions, 
including objects at high redshift. Thus, in Figure 14 we plot the model extinction resulting from 
several of our distributions over an extended wavelength range. 

In Figure 15, we plot the albedo and asymmetry parameter g = (cos 6) (i.e. the average value 
of cos 6*, where 9 is the angle through which radiation is scattered by dust) resulting from several 
of our model size distributions. 

Since Li & Draine (2001) find that the IR emission from dust in the diffuse ISM is best fit when 
6c ~ 6 X 10~^, we adopt this value for the Ry = 3.1 curves in Figures 14 and 15. For such a large 
6c, the 2175 A hump is almost entirely due to the very small carbonaceous grain population. If this 
is the case for the diffuse ISM, then it seems plausible that it also holds in denser regions; i.e., the 
decrease in the strength of the 2175 A feature with Ry might result entirely from the depletion of 
very small carbonaceous grains. Thus, we have also adopted the large-6c distributions for Ry = 4.0 
and 5.5 in Figures 14 and 15. 

3.1.3. Dust Along the Line of Sight to HD 210121 

Although the variation of the extinction curve with interstellar environment is fairly well 
characterized by the CCM parameterization, there are lines of sight for which the extinction deviates 
substantially from CCM. As a further test of the bare carbonaceous/silicate dust model, it is 
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important to seek size distributions which can reproduce the extinction along such sighthnes. The 
extinction observed toward HD 210121 (a sighthne passing through a high-latitude diffuse molecular 
cloud) has (1) an extremely small value of Ry = 2.1, (2) a 2175 A feature weaker than predicted by 
the CCM parameterization, and (3) a stronger-than-expected far-UV rise (see Figure 1 in Larson et 
al. 2000). This sightline therefore provides an opportunity to test the carbonaceous/silicate model 
and the functional forms used for our size distributions. 

Larson et al. (2000) used the maximum entropy method to construct size distributions for 
the grains toward HD 210121. We seek to reproduce the extinction toward HD 210121 (Larson 
et al. 2000; Larson, Whittet, &: Hough 1996; Welty Sz Fowler 1992) with size distributions of our 
simple functional form. We adopt the normalization given by Larson et al. (2000): Ay/Nn = 
3.6 X 10~^^cm^. In fitting the extinction, we adopt 100 points equally spaced in A^^ rather than 
in In A. We have found that this yields a better fit to the 2175 A hump and far-UV rise without 
compromising the fit quality in the infrared. Distribution parameter values are given in Table 2 
and the distributions and extinction fits are plotted in Figures 16 and 17, respectively. 

We are able to obtain acceptable fits to the extinction toward HD 210121 with values of be 
ranging up to 4 x 10~^, and reasonable size distributions for the carbonaceous and silicate grain 
populations. Our grain model successfully accomodates this line of sight with its extremely small 
value of Rv and deviation from the CCM parameterization. 

3.2. Dust in the Magellanic Clouds 

The metallicities in the Magellanic Clouds are substantially lower than in the Milky Way, and 
measured extinction curves toward stars in the LMC and SMC differ from typical extinction curves 
in the Milky Way. The LMC and SMC therefore offer opportunities to test the applicability of our 
grain model to low-metallicity cxtragalactic environments.^ 

Clayton et al. (2000) used the maximum entropy method to find graphite/silicate size distri- 
butions that accurately reproduce the extinction along various Magellanic Cloud sightlines. Here, 
we seek distributions of our simple functional form that reproduce the average extinction in the 
LMC (Misselt, Clayton, & Gordon 1999), the extinction in the LMC 2 area (Misselt et al. 1999), 

and the extinction in the SMC bar, along the line of sight to the star AzV398 (Gordon & Clayton 
1998). For A^-*^ < 3;um~^, the extinction is determined at only a small number of wavelengths. 
Thus, for the Magellanic Clouds, we evaluate the extinction at 100 wavelengths spaced equally in 
A"^, rather than in In A. 

The extinction normalization and elemental abundances are even more uncertain for the Mag- 
ehanic clouds than for the Milky Way. For the LMC, Koorneef (1982) found N{m)/E{B - V) = 



See Pei (1992) for an early extension of the MRN model to the Magellanic Clouds. 
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2.0 X lO^^cm-^ and Fitzpatrick (1985) found A^(HI)/£;(B - V) = 2.4 x lO^^cm-^. Averaging 
these results and taking Ry = 2.6 (the average for the 10 measured Ry values in Misselt ct al.'s 
sample), we adopt A{V)/Nu = 1.2 x lO'^^ ^^2^ -por the SMC, Martin, Maurice, & Lequeux (1989) 
found Nu/E{B - V) = 4.6 x 10^2 cm'^; with Ry = 2.87 (Gordon k Clayton 1998) this yields 
A{V)/N}i = 6.2 X lO^^^cm^. We take the abundance/depletion-limited values of Viot,g and Vtot,s 
to be reduced from their values in the Milky Way by a factor 1.6 for the LMC and 4.0 for the SMC 
(Gordon & Clayton 1998). 

Distribution parameters for which the extinction is best fit are given in Table 3. We also tabu- 
late the total grain volumes, normalized to the limiting values estimated in the previous paragraph; 
note that all of the LMC distributions use less than the estimated available amount of C and Si. 
Size distributions, extinction fits, and related quantities are plotted in Figures 18 through 23. 

Note the absence of the 2175 A feature in the SMC bar extinction curve (Figure 21), which 
implies the absence of very small carbonaceous grains. Recently, Reach et al. (2000) have detected 
PAH emission features in a quiescent molecular cloud in the SMC. Reach et al. point out that SMC 
extinction curve measurements are biased towards hot, luminous stars, so that very small grains 
may have been destroyed along these sightlines. 



Note from Table 1 that, in the Milky Way, the silicate volumes generally exceed the abundance/depletion- 
limited value, by « 10% when Ry = 5.5 to « 30% when Ry = 3.1, and the carbonaceous grain 
volume exceeds its abundance/depletion-limited value by 10% when Ry = 3.1. We would 
expect non-spherical grains to produce more extinction per unit grain volume than spheres, so 
that our violation of abundance constraints might be an artifact due to the use of only spheri- 
cal grains in our modelling. However, we have used the discrete dipole approximation (Draine & 
Flatau 1994; Draine 2000) to calculate extinction efficiencies for silicate grains of various shapes 
with a > 0.01/Lim and have found that the integrated extinction per grain volume, j (Cf,xt/y)dX 
integrated over A"''- G [0.35, 8.0]/im~^, varies only slightly with shape. 

Kim et al. (1994) sought to maximize the efficient use of grain volume by allowing more 
complicated size distributions. Although such an approach could lower the total amount of grain 
volume that we need to reproduce the observed extinction, we find such fine-tuning unappealing. 
It seems to us unlikely that nature has produced size distributions fine-tuned to maximize the 
extinction per volume over just the wavelengths where we are able to measure the extinction. We 
think it more likely that either the true elemental abundances in the ISM really are somewhat higher 
than in the Sun, or that the bare graphite/silicate model is inadequate in some more fundamental 
way. 



4. 



Discussion 



4.1. 



Abundances and Grain Models 
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Other well-developed models include composite, fluffy grains (Mathis 1996, 1998) and grains 
consisting of silicate cores covered by organic refractory mantles (Li &: Greenberg 1997). The recent 
discovery that the S.Afim aliphatic C-H stretch absorption feature toward Sgr A IRS7 is unpolarized 
(whereas the 9.7//m silicate absorption feature toward Sgr A IRS3 is polarized) may rule out the 
core-mantle model (Adamson et al. 1999), although model calculations of the relative polarization 
in these features have not yet been carried out for the core-mantle model, and the silicate feature 
polarization has yet to be measured for IRS 3 itself. 

Mathis (1996) found that a mixture of composite grains (consisting of small silicate and amor- 
phous carbon grains and ^ 45% vacuum), small graphite grains, and some small silicate grains 
could reproduce the observed extinction while incorporating C, Si, Fe, and Mg with substantially 
sub-solar abundances. However, there are some difficulties with this model. First, Mathis adopts 
dielectric functions for the composite grains using effective medium theory, calculates extinction 
cross sections for spheres, and then multiplies the cross sections by a factor 1.09, to account for en- 
hancements in extinction due to non-spherical shapes. The flnal step must be viewed with suspicion, 
since it fails for compact silicate grains. 

Also, Mathis used the optical properties of "Be" amorphous carbon from Rouleau & Mar- 
tin (1991). Schnaiter et al. (1998) have pointed out that the derived optical properties, while 
possibly correct, are unproven, since the adopted description of the sample geometry as a con- 
tinuous distribution of ellipsoids is so simplistic that substantial errors can result. Furthermore, 
"Be" amorphous carbon is much more absorbing at long wavelengths than various forms of hydro- 
genated amorphous carbon, and this absorption provides most of the extinction for < 3/Ltm~^ 
in the Mathis (1996) composite model. ^ Furton, Laiho, &; Witt (1999) have performed laboratory 
studies of hydrogenated amorphous carbon, and flnd that such grains can reproduce the observed 
3.4^m absorption feature if the degree of hydrogenation is rather large (~ 0.5 H/C). There is very 
little visible/IR continuum absorption in this case. Thus, the composite model docs not simul- 
taneously provide enough long-wavelength extinction and 3.4//m absorption. Of course, the bare 
graphite/silicate model does not account for the 3.4//m absorption either. ■'^'^ 

Although the bare graphite/silicate model apparently requires higher abundances of C, Si, Fe, 
and Mg than are generally thought to be available in the ISM, it would be premature to abandon it. 
The true interstellar abundances are not yet known, and the alternatives have difficulties too. Fur- 
ther progress in dust modelling will require the determination of dielectric functions for amorphous 
carbons with a range of degrees of hydrogenation, over the full range G [0.35, 8.0]/im~^, as well 
as detailed modelling of how the extinction per unit volume varies depending on grain geometry. 



^Dwek (1997) has argued that the fluffy grain model employing "Be" amorphous carbon produces too much IR 
emission compared with the COBE data (Dwek et al. 1997). 

^"To accomodate the 3.4/um feature, the graphite/silicate model must be extended to include aliphatic hydrocar- 
bons, possibly within hydrogenated carbon coatings on the large graphitic grains. 
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4.2. Observed Size Distribution of Interstellar Grains Streaming Through the Solar 

System 

Recently, Frisch et al. (1999) have presented a grain mass distribution for the local interstellar 
medium (LISM) , derived from the measured rate of impact of interstellar grains with detectors on 
the Ulysses and Galileo spacecraft; we reproduce their data points in Figure 24. We also show 
mass distributions as derived here from fitting extinction, for (iiy, 10^ 6c) = (3.1, 3.0), (4.0, 2.0), 
and (5.5, 1.0). We adopt nn = 0.3 cm~^, as recommended by Frisch et al. Note that none of our 
distributions resemble the Frisch et al. result. The steep drop in the Frisch et al. distribution at 
small masses probably reflects the exclusion of small grains from the solar system - smaller grains 
are more tightly coupled to the magnetic field and are less likely to penetrate the heliosphere to 
within ~ 5AU of the Sun (Lindc & Gombosi 2000). However, the large amount of mass in large 
grains in the Frisch ct al. distribution is hard to fathom. The error bars on the Frisch et al. data 
(not shown in Figure 24) are large; further observations of interstellar dust entering the solar system 
would be of great value. 

If the Frisch et al. result is confirmed, then there are two possibilities. If the region through 
which the solar system is now passing contains a truly representative dust-gas mixture, then a 
dramatically different grain model would be required. It is difficult to envision a grain model which 
could simultaneously account for the interstellar extinction law, be consistent with interstellar 
elemental abundances, and reproduce the Frisch et al. size distribution. Alternatively, it could be 
the case that size-sorting and gas-grain separation occur on small scales in the ISM, and that the 
region through which the solar system is now moving happens to have an unusual concentration of 
large grains. 

4.3. Conclusions 

The simplest interstellar dust model consists of a population of carbonaceous grains and a 
separate population of silicate grains. In the original development of this model by MRN, the 
grain size distribution was chosen so as to reproduce the observed extinction for lines of sight with 
Rv ~ 3.1. The observation of relatively short- wavelength infrared emission from dust implies that 
there are substantial numbers of very small (mainly carbonaceous) grains, smaller than the lower 
cutoff size of the MRN distribution. Furthermore, the extinction curve has been found to vary 
substantially depending on the interstellar environment through which the starlight passes; thus, 
there is no single grain size distribution which applies in all environments. By finding carbona- 
ceous/silicate grain size distributions which contain sufficient very small grains to account for the 
observed infrared emission (Li & Draine 2001), and which reproduce the observed extinction for a 
wide range of environments, we have demonstrated that the simplest dust model remains viable. 

Although difficulties remain, they are no more severe than the difficulties with other, more 
complicated, models. These difficulties include the requirement of somewhat super-solar abun- 
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dances of the dust constituent elements, the lack of a 3.4/xm absorption feature in a model in which 
all of the C is in graphite or PAHs, and the gross disparity between the derived grain size distri- 
butions and that inferred by Frisch et al. (1999) for dust in the local ISM. Additionally, there is 
evidence from depletion patterns that metallic Fe or Fe oxides are an important dust component 
(Sofia et al. 1994; Howk et al. 1999). The observed 90 GHz emission from interstellar dust appears 
to rule out a substantial metallic Fe component (Draine &; Lazarian 1999), but oxides such as FeO 
or magnetite Fe304 are not excluded. Dielectric functions for candidate Fe oxides are needed to 
investigate such grain models. 

Finally, the variation in the grain size distribution with environment seems to indicate that 
small grains coagulate onto large grains in relatively dense environments, as expected (Draine 1985; 
Draine 1990). Presumably, mass is returned from large to small grains via shattering during grain- 
grain collisions in shock waves. (Mass is also returned to the gas via sputtering processes.) Wein- 
gartner &; Draine (1999) found that the observed elemental depletions in the interstellar medium 
could be due to accretion onto grains if the timescales for matter to cycle between interstellar phases 
are ~ 10^ yr. It remains a mystery how two separate grain populations - carbonaceous grains and 
silicate grains - could remain distinct after evolving through many cycles of coagulation, shattering, 
accretion, and erosion; perhaps they do not. 

While real grains are undoubtedly more complex, the graphite/silicate model for dust in dif- 
fuse clouds is clearly-defined, and consistent with observations of interstellar extinction in the Milky 
Way, LMC, and SMC (as demonstrated in the present work) and infrared emission (Li &i Draine 
2001). While the model does not explicitly account for the 3.4/xm feature or the relatively weak 
diffuse interstellar bands (Herbig 1995), these could conceivably be accomodated by modest mod- 
ifications of or extensions to the basic graphite/silicate model. The "extended red emission" from 
interstellar dust (Witt &; Boroson 1990) could also perhaps be due to a minor modification of 
the basic graphite/silicate model (e.g., a hydrogenated amorphous carbon coating; Witt Sz Furton 
1995). 

Until a more compelling grain model is available, we recommend the use of the simplest one, 
specified by the size distributions found here and optical properties given by Draine &: Lee (1984), 
Laor & Draine (1993), and Li Sz Draine (2001). In particular, we favor the distributions with 
relatively large 6c (Li Sz Draine 2001), for which the very small carbonaceous grain population 
entirely accounts for the 2175 A hump in the extinction curve. 

This research was supported in part by NSF grant AST-9619429 and by NSF Graduate and 
International Research Fellowships to JCW. We are grateful to Eli Dwek, Aigen Li, and John 
Mathis for helpful discussions and to R. H. Lupton for the availability of the SM plotting package. 
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Note added in proof — Li Sz Draine (2000) have recently found that the nondetection of the 

10/7,m silicate feature in emission from diffuse clouds docs not strongly constrain the ultrasmall 
silicate grain population, since the W/im feature may be hidden by the dominant PAH features. Li 
&; Draine estimate that as much as ~ 20% of the interstellar Si could be in grains with a < 15 A. 



This preprint was prepared with the AAS IM^^X macros v5.0. 
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Fig. 1. — The error function versus the siUcate cutoff parameter, ac,i 
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a (/am) 

Fig. 2. — Case A grain size distributions for Ry = 3.1. The values of be are indicated. The heavy, 
soHd Hues are the MRN distribution, for comparison. Our favored distribution has 6c = 6 x 10~^ 
(see text). 



2008.2.1.1450: DRAFT 
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Fig. 3. — Same as Figure 2, but for Ry = 4.0. Our favored distribution has 6c = 4 x 10~^ 
text). 
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Fig. 4. — Same as Figure 2, but for Ry = 5.5. Our favored distribution has 6c = 3 x 10 
text). 
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Fig. 5.— 



Case B size distributions for Ry = 4.0. 
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Fig. 6.— 



Case B size distributions for Ry = 5.5. 
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1 

A-i (/xm-i) 

Fig. 7. — The average "observed" extinction AqBs and the extinction resulting from our case A 
models for {Ry, lO^fec) = (3.1, 6.0), (4.0, 4.0), and (5.5, 3.0). The curves for Ry = 4.0 (5.5) are 
scaled down by a factor lO'''^ (10°'^), for clarity. 
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2 4 6 8 

A-i (/xm-O 

Fig. 8. — The extinction curve ^mod resulting from the grain distribution of equations (4) and 
(5), with parameters optimized to fit ^obs (sec text) for Ry = 3.1 (also shown), for be = 0.0 and 
6.0 X 10~^. The contributions from the three grain distribution components are also shown. 
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2 4 6 8 

A-i (Mm-i) 

Fig. 11. — Same as Figure 8, but for Ry = 4.0, 6c = 0.0 and 4.0 x 10~^, and fixed total grain 
volumes T4ot,g = 2.3 x lO-^^cm^H-^ and Ftot.s = 3.9 x lO-^^ cm^ H-^ 
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Fig. 13. — The extinction fit error function Xi (§2-5) as a function of be, the C abundance in the 
log-normal grain population, for three values of Ry- 
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Fig. 15. — Albedo and asymmetry parameter g = {cos 6) for various size distributions. 
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Fig. 16. — Grain size distributions for HD 210121. 
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2 4 6 8 



Fig. 17.- Same as Figure 8, but for the extinction along the hne of sight to HD 210121 and 
6c = 0.0 and 4.0 x 10~^. Note the difference in vertical scale from Figure 8. 
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Fig. 18. — Grain size distributions for the LMC. The values of 6c are indicated; "A" denotes 
distributions constructed to fit the average extinction in the LMC and "2" denotes distributions for 
the LMC 2 area. 
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Fig. 19. — Same as Figure 8, but for the average extinction for the LMC and be = 0.0 and 2.0 x 
Note the difference in vertical scale from Figure 8. 
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Fig. 20. Same as Figure 8, but for the LMC 2 area and be = 0.0 and 1.0 x 10"^. Note the 
difference in vertical scale from Figure 8. 
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Fig. 21. — Upper panel: Size distribution for the SMC bar, with be = 0.0. Lower panel: The 
corresponding extinction fit; curve types are the same as in Figure 8. 
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Fig. 22. — Model extinction curves extended to short wavelengths, for Magellanic Cloud environ- 
ments. 
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Fig. 23. — Albedo and asymmetry parameter g = (cos 9) for Magellanic Cloud environments. 
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Fig. 24. — Mass distribution for grains in the local ISM determined by Prisch et al. (1999) (tri- 
angles). Mass distributions for the size distributions of §3 are also shown; the sharp drop at 
m ~ 3 X 10~^^ g corresponds to the rapid drop in silicate grain abundance at a ~ 0.3/xm. 
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X 


10- 


14 


0.933 


1.145 


0.037 


0.148 


0.060 


5.5 


0.0 


A 


-2 


35 


-0.668 


0.148 


1.96 


4.82 


X 


10" 


14 


-1.57 


1.10 


0.198 


4.24 


X 


10- 


14 


0.889 


1.076 


0.034 


0.110 


0.043 


5.5 


1.0 


A 


-2 


12 


-0.670 


0.0686 


1.35 


3.65 


X 


10- 


13 


-1.57 


1.25 


0.197 


4.00 


X 


10- 


14 


0.848 


1.078 


0.034 


0.115 


0.043 


5.5 


2.0 


A 


-1 


94 


-0.853 


0.0786 


0.921 


2.57 


X 


10- 


13 


-1.55 


1.33 


0.195 


4.05 


X 


10- 


14 


0.804 


1.095 


0.032 


0.118 


0.044 


5.5 


3.0 


A 


-1 


61 


-0.722 


0.0418 


0.720 


7.58 


X 


10- 


13 


-1.59 


2.12 


0.193 


3.20 


X 


10- 


14 


0.768 


1.118 


0.033 


0.128 


0.049 


4.0 


0.0 


B 


-2 


62 


-0.0144 


0.0187 


5.74 


6.46 


X 


10- 


12 


-2.01 


0.894 


0.198 


4.95 


X 


10- 


14 






0.011 


0.042 




4.0 


1.0 


B 


-2 


52 


-0.0541 


0.0366 


6.65 


1.08 


X 


10- 


12 


-2.11 


1.58 


0.197 


3.69 


X 


10- 


14 






0.011 


0.043 




4.0 


2.0 


B 


-2 


36 


-0.0957 


0.0305 


6.44 


1.62 


X 


10- 


12 


-2.05 


1.19 


0.197 


4.37 


X 


10- 


14 






0.011 


0.042 




4.0 


3.0 


B 


-2 


09 


-0.193 


0.0199 


4.60 


4.21 


X 


10- 


12 


-2.10 


1.64 


0.198 


3.63 


X 


10- 


14 






0.011 


0.044 




4.0 


4.0 


B 


-1 


96 


-0.813 


0.0693 


3.48 


2.95 


X 


10- 


13 


-2.11 


2.10 


0.198 


3.13 


X 


10- 


14 






0.017 


0.056 




5.5 


0.0 


B 


-2 


80 


0.0356 


0.0203 


3.43 


2.74 


X 


10- 


12 


-1.09 


-0.370 


0.218 


1.17 


X 


10- 


13 






0.017 


0.092 




5.5 


1.0 


B 


-2 


67 


0.0129 


0.0134 


3.44 


7.25 


X 


10- 


12 


-1.14 


-0.195 


0.216 


1.05 


X 


10- 


13 






0.017 


0.088 




5.5 


2.0 


B 


-2 


45 


-0.00132 


0.0275 


5.14 


8.79 


X 


10- 


13 


-1.08 


-0.336 


0.216 


1.17 


X 


10- 


13 






0.017 


0.085 




5.5 


3.0 


B 


-1 


90 


-0.0517 


0.0120 


7.28 


2.86 


X 


10- 


12 


-1.13 


-0.109 


0.211 


1.04 


X 


10- 


13 






0.017 


0.082 





''See equations (4) and (5). In all cases, we take ac,B = 0.1/um. 
^Rv = A{y)/EB-v, ratio of visual extinction to reddening 

"^C abundance in double log-normal very small grain population (see equations 2 and 3) 

"^Total grain volumes in the carbonaceous and silicate populations, normalized to their abundance/depletion-limited values (2.07 x 10"^^ and 2.98 x 10"^^cm^H 
respectively) 

®Xi = SiG'^^obs ~ In Aniod)^/"'! I points equally spaced in InA 

^Xl = Ei(ln^obs - InAniod)^ 

=xl+ 0.4(Vg - 1)1-5 -t- 0.4(^3 - l)i-« 



Table 2. Grain Size Distribution Parameter Values for HD 210121^ 



106 be'' 




/3g 






Cg 




as 






Cs 






^4-= 






X^' 








(/im) 


(/im) 










(/im) 
















0.0 


-2.22 


-0.0960 


0.00544 


0.651 


1.71 X 10" 


-10 


-1.96 


-5.23 


0.0999 


2.32 X 10" 


■ 12 


0.752 


1.407 


0.071 


0.080 


0.175 


1.0 


-2.18 


-0.0818 


0.00551 


0.614 


1.28 X 10" 


-10 


-1.98 


-5.25 


0.105 


1.99 X 10" 


■ 12 


0.745 


1.415 


0.070 


0.078 


0.177 


2.0 


-2.04 


-0.137 


0.00731 


0.566 


5.37 X 10" 


-11 


-1.96 


-6.05 


0.110 


1.97 X 10" 


-12 


0.736 


1.423 


0.069 


0.077 


0.179 


3.0 


-1.87 


-0.190 


0.00911 


0.492 


2.40 X 10" 


-11 


-1.94 


-6.99 


0.112 


2.09 X 10" 


■ 12 


0.726 


1.428 


0.072 


0.082 


0.184 


4.0 


-1.69 


-0.264 


0.0126 


0.449 


8.60 X 10" 


-12 


-1.90 


-9.22 


0.119 


2.26 X 10" 


■ 12 


0.715 


1.442 


0.077 


0.088 


0.194 



^Scc equations (4) and (5). In all cases, we take Oc.s = 0.1/im. 

'^C abundance in double log- normal very small grain population (see equations 2 and 3) 

"^Total grain volumes in the carbonaceous and silicate populations, normalized to their abundance/depletion-limited values (2.07 x lO"^'^ and 2.98 x 

10"'^^ cm^ , respectively) 

'^Xi = Z]i(lii^obs — 111 ^mod)^/"'! I for 100 points equally spaced in 

"xi = Ei(iii^obs-iiiA^od)^ 

^X^ = xi + 0.4(14 - 1)1-5 + OAiV, - 1)1-5 



Table 3. Size Distribution Parameter Values for the Magellanic Clouds'^ 



Environment 




"g 


/3g 


(^m) 


ac,g 
(^m) 






as 




a-t,B 
(/im) 


Cs 






v." 




xV 




LMC avg 


0.0 


-2.91 


0.895 


0.578 


1.21 


7.12 X 10" 


-17 


-2.45 


0.125 


0.191 


1.84 X 10" 


-14 


0.401 


0.675 


0.025 


0.069 


0.025 


LMC avg 


1.0 


-2.99 


2.46 


0.0980 


0.641 


3.51 X 10-^ 


-15 


-2.49 


0.345 


0.184 


1.78 X 10- 


-14 


0.330 


0.687 


0.018 


0.033 


0.018 


LMC avg 


2.0 


4.43 


0.0 


0.00322 


0.285 


9.57 X 10" 


-24 


-2.70 


2.18 


0.198 


7.29 X 10" 


-15 


0.279 


0.758 


0.016 


0.019 


0.016 


LMC 2 


0.0 


-2.94 


5.22 


0.373 


0.349 


9.92 X 10" 


-17 


-2.34 


-0.243 


0.184 


3.18 X 10" 


-14 


0.263 


0.753 


0.025 


0.043 


0.025 


LMC 2 


0.5 


-2.82 


9.01 


0.392 


0.269 


6.20 X 10" 


-17 


-2.36 


-0.113 


0.182 


3.03 X 10" 


-14 


0.252 


0.765 


0.022 


0.037 


0.022 


LMC 2 


1.0 


4.16 


0.0 


0.342 


0.0493 


3.05 X 10- 


-15 


-2.44 


0.254 


0.188 


2.24 X 10" 


-14 


0.206 


0.820 


0.012 


0.014 


0.012 


SMC bar 


0.0 


-2.79 


1.12 


0.0190 


0.522 


8.36 X 10- 


-14 


-2.26 


-3.46 


0.216 


3.16 X 10" 


-14 


0.254 


1.308 


0.017 


0.019 


0.027 



^Scc equations (4) and (5). In all cases, we take ac,s = 0.1/im. 

''C abundance in double log-normal very small grain population (see equations 2 and 3) 

''Total grain volumes in the carbonaceous and silicate populations, normalized to their abundance/depletion-limited values (1.29, 1.86, 0.518, and 0.745 x 10~^^ cm^ H~ 
for carbonaceous in LMC, silicate in LMC, carbonaceous in SMC, and silicate in SMC, respectively) 

''Xi = X^iC'i^obs — In Aniod)^/ff 1 for 100 points equally spaced in X~^. 

''xi = Ei(ln Aobs - In A^od)^ 

^X^ = Xi + 0.4(^8 - 1)1-5 + 0.4(14 - l)l-5 



